THE DISTRIBUTION OF RANDOM EVOLUTION IN ERLANG 
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j Abstract. In this paper we study a one-dimensional random motion by hav- 

jrt ■ m g a general Erlang distribution for the sojourn times of the switching process 

and we obtain solution of the four order hyperbolic PDE for 2-Erlang case. 



On 



1. Introduction 



(r , • In the paper [1] we studied a one-dimensional random motion with the m-Erlang 

P I | distribution between consequent epochs of velocity alternations. Let f(t,x) be the 

(■h ■ probability density function (pdf) of a particle position at time t, provided that it 

-^ exists. We obtained the following higher order hyperbolic equations for f(t,x) 

m -v^ + x) (^ + ^ + a) fit,,)- A-/ (i ,x) = 0, (1) 

where v > is the velocity of the particle and A is the parameter of the m-Erlang 
distribution. It is assumed that a particle started at x — and hence, f(0,x) = 

The pdf f(t,x) can be represented in the following form f(t,x) = f c (t,x) + 
in | fs(t,x), where f c (t,x) is the absolute continuous part and f s (t,x) is the singular 

part w.r.t. Lebesgue measure on the line. 

Lemma 1. TTie singular part f s (t, x) of the pdf f(t, x) is of the following form 

O m-l 

/,(*, x/v) = S(t - x/v)e~ xt J2 (A*)V*1. (2) 

i=0 



k> ' Proof. It is evident that for £ = x/v the pdf f(t,x) has the singularity given 

5_i , by Eq.(2). Let us show that for t > \x/v\ the pdf f(t,x) has no singularity w.r.t. 

Lebesgue measure on R. Denote by Vk the random event "fc velocity alternations 
occurred". For Aa; — [x,x + A], A > 0, consider 

Pp (x{t) g Ax) =5^P(a:(t) G Aar.i^), 
fe>i 

which is the probability of the event that at least one alternation occurred and 
x(t) G Ax. Let us show that for each t > there exists a constant Ct < oo such 
that 

wmM <Ct . 

x Ax 
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Denote by 9k, k > 1, time between (k— l)th and fcth velocity alternations. Recall 
that 9k, k > 1 arc independent m-Erlang distributed random variables. It is easily 
verified that 

(k / k \ k \ 

53(-l)' +1 0^ + (-l) fe t - 2^v U Ax, ^ ft < i 
j=l \ i=l / i=l / 

53(-l) i+1 0» - (-l) fe jS* ) V G As - (-l) k vt, J2 9 % < t 

\i=l i=l / i=l J 

= ^ P 1 2v(9 1 + 3 + . . . + 9 2l+l ) e Ax - vt, ]T $i < t J 

Z>0 V »=1 / 

= J2 P [ -2v(9 2 + 04 + . . . + 2;+2 ) € Ace + vt, 53 ft < * J 

l>0 V »=1 / 

< sup 53 p ( 2w 5Z ft*- 1 e Aa; ' 2 « 5Z ^ 2i < wf _ x ) 
+ sup 53 p ( -2v 53 ft* e Aa: ' 2w X! ^ 2i+i < ^^ + ^ ) • 

2 ;>o V i=l »=1 / 

Since | a; | < vt and for every m > 1 the pdf p m (x, A) of the m-Erlang distribution 
with the parameter A satisfies p m {x, A) < A, we have 

53 P(2v(6i + 6> 3 + • • • + 02i-i) G Ax, 2v(9 2 + 4 + • • • + fti) < vt - ») 



/>i 



<^Y / P(92 + 9 i + ... + 9 2l <t) (3) 



2w 
;>i 



Since 0^ is m-Erlang distributed we have for 2lm + 1 > t 



llm 



P(0 2 + 04 + . . . + 92i < t) < e At - 53 ^- e- At < 



(At)' \ , t ^. (Atf^e 



2lm+l —Xt 



i\ I ~ 2lm\(2lm + I - At)' 

Therefore, taking into account (3), there exists a constant A t such that 

sup 53 p ( 2v53#2i-i e Ax, 2^53^ < vt-x ] < A t A. 
11 ;>i V »=i »=i / 

In much the same way, we can show that there exists a constant B t such that 

CI 2Z-1 \ 

-2t> 53 hi e Ax ' 2u Y foi < vt + x) < B t A. 
j=l i=l / 

Putting Ct ~ A t + B t , we conclude the proof. 

Corollary 1. The absolute continuous part / c (t, x) of the pdf f(t, x) satisfies Eq. (1) 
/or* < If I- 
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Now let us study the behavior of the continuous part f c (t, x) close to lines t = ±-. 
Lemma 2. For m > 2, we have 

,. P{0 < t - x(t) < el X m t m - 1 e- Xt 

hni = ; rr- , 

do e 2(m-l)! 

lim P{t + ^ )<£} = 0. 



£-i-0 

Proof. It is easily verified that 



pi 
P{0<t-x(t) <e} = PU-^<6 1 <t\+ P{e 3 >t-u,6 2 < E -M Edu\+o(s), 

where 9i,i — 1,2,3 are independent rn-Erlang distributed random variables with 
the parameter A. Since L P($3 > t — 11,82 _• §,#1 € du) = o(e), passing to the 



limit, we get 






,. P{0<t-x(t) <e} 






siO € 






_1 ™ s \{h i - J 


-e A f 


V t=0 


X m t m-l e -Xt 







2(m-l)! ' 

Similarly, P {£ + x (i) < e} = -P {# 2 > t — §, #1 < §} + o (e) and as it easily seen 
that 

lim P{t + x(t)<e} = Q 

eiO £ 

The case where m — 1 will be considered below as an example. 
We will seek solutions of Eq.(l) among functions which continuous part f c (t,x) 
satisfies the following conditions 

lim /„(*,*) -lim P{0< ^ (t)<£} , lim/ cM lim P ^ + * (t) < e} _ 

xft e|0 £ xl—t £j,0 £ 

By applying the transformation / (t, x) — e xt g (t, x) and changing the variable 
y = -, we reduce Eq.(l) to 

W'W) 9c(t,y)-X 2m g c (t,y) = 0, (5) 

with the singular part g s (t, y) = (j]™ q 1 ^~) 5 ( l ~ v) ■ 

In the sequel we assume, without restricting the generality, that A = 1. By 
introducing the function f (t,y,z) — e z g c (t,y), we reduce Eq.(5) to the following 
equation 

d 2 d 2 \ m d 2m 

We will seek solutions of this equation by using theory of differentiable functions 
on commutative algebras [2]. 
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2. Main results 



Let Aq be an 2m- dimensional commutative algebra over R, assume that the set 
eo, ei, . . . , e2 m _i is a basis of ^4o with the Cayley table: 



where i (B j = i + j {mod 2m) . 

Algebra Ao has the following matrix representation: 

e k -> P fe = P*, 

where Pi = [Pij] 2mX 2m' P«+i = 1 for < i < 2m - 1, p 2m o = 1 and p^- = for the 
rest of i, j. 
Let us put 

r = ei, I = 0,1, ... , 2m — 1, 

r{=e;isins, I = 0, 1,. . .,2m— 1, 

t 2 = e; cos s, I = 0,1, ... , 2m — 1, 

T 2fe = e i cos ^ s i T 2fe+i = e '* sin(/c + l)s, I = 0,1, ... , 2m — 1, 

k = 0,l,2,.... 

It is easily that T 2 °„T 2 ° fe = § (r 2 ° (ri _ fe) + T 2 ° (n+fe) ) , n > fc, 

Ji J2 _ 1 fji®h Ji®k 



'2n+l'2k+l — 2 \ 2(n-k) 2(n+k)J ' 

T h T h_l( T h®l2 , Ii©Ja \ >r. 

T 2n+l T 2fc _ 2 \ T 2(n-fe) + l + T 2(n+fc) + lJ ' ft - ft ' 



Let us introduce the following algebra 

( +oo 2m-l 
-"■ = \ 2^i Z^i \ a 2k T 2k + a 2fc+l r 2fc+lj 



,fe=0 1=0 
^+co Y^2rn — 1 /l / |2 i ^ |2 



where LtoLtO (J a 2fc| + Kfc+l| J < +°°- 

It is easily verified that A is commutative. 

We consider the subspace B = {ciqtI + cl\t\ + a 2 TQ | a, € R} of the algebra A. 
Let us introduce the function f : B — » A (f (£, y, z) = / (ei (teoss + yisins ) + z)) 
as follows 

+oo 2m- 1 

f{t,y,z) = Y^ Yl ( v 2k(^y^)4k + v 2k+i(^y^) T 2k+i)- 

k=0 1=0 

The function f is called B/A diffcrentiable at xo € B if there exists f (xo) € A 
such that for any h <G B 

f (x )h=lim f(X0 + £h) ' f(X0) 
In [2] proved that if f is P/A diffcrentiable, then 

— f = eiCOSs — f (7) 

at oz 

and 
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— f = eiisins — f. (8) 

dy oz 

In this case all v l 2k (t, y, z) are solutions of Eq.(6). Indeed, 

In the sequel we denote by e the element e\. 

We will seek a solution of Eq.(5) in the following form 

g c (e (teoss + yisins )) = e e (* coss +^ sins ) 

Since / (e (teos s + yism s ) + z) = g c (e (tcos s + yisin s ))e z we have 

v{ (t, y,z) = u{ (t, y) e z ,1 = 0, 1, . . . , 2m - 1, k = 0, 1, 2, . . . , 

where g c (t, y) = £+^ EST* { u \k (*. v) T L + u \k+i (*> v) T 2fc+i) • 

Therefore, we obtain functions u l (t,y) for t > |j/| from the following equation 

2m- 1 2m- 1 

^ «o(t,y)7^= 51 «o(*.I/)e' 



/-() /-I) 

?y 2 



= — f e e(tcoss +yisins ) ds = J (eiy/y 2 - tA = I (eV^ - 



where /& (resp. Jfc) is the modified Bessel (resp. Bessel) function of the first kind 
and kth order [4]. 

where Ik (resp. Jk) is the modified Bessel (resp. Bessel) function of the first 
kind and fcth order [4]. 

It follows from Eqs.(7),(8) the following Cauchy-Riemann type conditions 

9 (0i ! i 
dt ° 2 2 ' 



d i®i 1 . 
dt 1 2 3 

— u 1 ® 1 -u l +-u 

dt 2 - U "+2 



,i 



~gj. u 2k-l — 2 V M 2fe-3 + u 2k+l) i 



Q ' \T7 1 -I 

-^- = 2( u 2(fe-i)+4(fe + i)) (9) 

fc = z, o, . . . ; 



and 



9 J©1 

9/° 


= - 


4«i. 


flu 1 © 1 

dy 


u l 


-^«i, 


dy 


= - 


~t4, 
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I 



ua 2k+l _ l ( I 

dy ~ 2 v 2fc U2 ( k + 2 '))> 
On 1 ® 1 I 



fiy 



w 



2fe-l u 2fc+3j i 



(10) 



fc=l,2,.... 

By using Eqs.(9),(10) and functions it (t,y), we can obtain recurrently function 
it' fe (t, y) for any fc > 1, which will be used for solution of Eq.(l). 

In the sequel, unless otherwise specified, the case where m = 2 is studied. In 
this case f s (t, y) is of the following form f a (t, x) = e~ l (1 + 1) 6 (t — x) and hence, 

g s (t,y) = (l + t)5(t-y). 

Algebra Aq is as follows 

Aq = {a + eib + e 2 c + e 3 d \ a, b, c, d E M. } , 

where the basis e; = e , I = 0, 1, 2, 3, and e has the following matrix representation 



P 



Therefore, we have Tq = 1, t^j. = cos fcs , r^ = e'cos fcs , T2 fe+1 = e'zsin (fc + 1) s , 
/ = 0,1,2,3, k = 0,1,2,... 

Taking into account that g c (e (tcoss + yisins )) = e e (' coss +v islns ), we have 






1 














1 














1 


1 












^(t,y)+eu 1 (t,y)+e 2 u 2 (t,y)+e 3 4(t,y) = — e^ tcoss +^ sins Us = I (e^^) . 



It is easily seen that 



I Uy/W=tf 



h (ye^f) + Io (t v^v) 2 / /o (v^^F) - lo (i Vt^~7) 



to ( V^^7) + Jo (y^^7) 2 h (V^?) - Jo (/?=¥) 



Therefore, for t > \y\, we have Uq (t, y) — Uq (t, y) — and 



«o (*,») = 



«o (*, y) = 



/o(v^" 


- y 2 ) + Jo (JW~- 


-y 2 ) 


2 


/o [Vt 2 


- v 2 ) - Jo [V?- 


-y 2 ) 



It follows from the first two equations of (10) that 
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1 r, 9 2 



h (V^¥) - Jo (V^^¥) 



dy 



V^^V 2 



ii V^?+/iK/^r 



v 3 - -2—7/° 

«i- 2 a /o 



/o (v 7 ^^) + ^o (v^^ 



dy 



V* 



h(V^ 2 ~y 2 )-Ji(Vt 2 ~^y : - 



Then it follows from the Cauchy-Riemann type conditions (9) that 



«2 (t, y) = 2 



ftig (t, g) g ( 7 ° (v^^) - j (v^^)) 



at 



a/ 



L=(/ 1 (v^^)+^ 1 (v^^)) 



dy 



«2 (*,») = 2- 



Similarly, for it 3 we have 



K (*,„) _ 5 [ 7 ° (v^F) + Jo (y/W=?) 



at at 



2iy 



m§ =2^-ui = 2- 

3 at x at 



(*2 - y 2)- 






+ ^T ( 7 ° (V^ 3 ^) + J 2 (v^V*) + Jo (v/t 1 ^) - J 2 (v 7 ^^)) 



<2-y 
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2ty 



ui =2— it? = 2— 

3 at 1 dt 



Vt^y 2 



h(V^^y 2 )-Ji(V^^y i 



(f -y 2 Y 



hlV^f )-Ji ^ 



+ 



2ty 

t 2 - ?/ 



(/o (v^^y 2 ) + h (yw^?) - Jo (yw^) + j 2 (yw^)) 



It is easily seen that u\ = wf = 0. 
Next, it follows from (9) that 



u° A =2 



du\ 



2w = 2 



a t 



-2y 2 



(/! (v^V*) + Ji (V^F)) - 2ug 



+^2 ( j o (v 7 ^^^) + ^ (v^^) + Jo (yp^¥) - h (y^^f)) 



ui=2 



~~dT 



-2y 2 



(h(V^^)-Ji(V^^)) 



\/(t 2 -y 2 ) 3 

+^2 ( 7 « (v^s) + h (v^s) Jo (v^^) + J2 (y^^f)) 

-h [yp^y 2 ) + Jo [yp^y 2 ) ■ 



Also it is easily verified that u\ = u\ = 0. 

By using well-known integrals for Bessel functions [3-5], we have 

/u dy = sinht + sint , / u dy = sinht — sint , / u x dy = / Uidy = 0, 
-t J-t J —t J—t 

ft ft Q u 2 / ft 

J u\dy=2j -^dy = 2(—J u\dy - u%{t,t) - u 2 (t,-t) 
=2cosht — 2cos£ , 

J uldy=2J fb£dv = 2\jtj u , dy-u° (t,t)-u , (t,-t) 

=2cosht + 2cos£ - 4. 

As an example, we obtain the pdf for the case, where m — 1. For this case ei = 1 
and hence, we can consider functions X^=o u k (A y)> k = 0,1,2, .. . as solutions of 
Eq.(5) for m = 1. 



THE DISTRIBUTION OF RANDOM EVOLUTION... 



For t < \y\ consider the function g (t, y) = g c (i, y) + g s (£, y) of the following form: 
\{u° Q {t,y) + ul(t,y)) 
/o(v^?) t 



9c (t, y) =2 («8 (*. v) + u l (*> y)) + 4 ( w i (*> v) + u l (*> y) + «2 (*, y) + «2 (t, y)) 



2^*^ 



y 



iliiVP^y* 



and g a (t,y) = 6(t-y). 

It is easily seen that function g c (t, y) is a solution of the equation for t < y 

d 2 d 2 \ 

0t2-Wj 9{t,y) ~ 9{t,y)=0, (U) 

In addition, we have lim^ g c (t, y) = | (1 + t) and lim^j-t g c (£, y) = 5 ■ 
To avoid cumbersome calculations we put v = 1. 
Therefore, / (£, x) = e~*g (£, x) is a solution of the equation: 

|-^ + 1 )(l + i + 1 ) /c(i ' X) - /c( ^ )= °' (12) 

/ s (i,x)=(5(i-x)e- t . 
In addition, / c (£, x) satisfies the following conditions: 

lim/ c (*, x) = - (e _t + ie - *) , Jim f c (t, x) = -e - *, 

and for all i > we have /_, f (t, x)dx = 1. 

For a small e > consider the probability P {0 < t — x (£) < e}. 

Let us verify that lim x1 - t / c (t, x) = lim^o P ~ x , i.e., 

,. P{0<t-x{t) <e} 1 , _ t ± _ n 
hm — — = - (e * + te ') . 

e|0 £ 2 V ' 

Indeed, it is easily seen that 

P {0 < t - x (t) < e} = P {* - £ - < $! < t} + / P {(9 3 > t - «, 62 < |, 0i g du} 

+o(e), 

where #j, i = 1,2,3 are independent exponentially distributed random variables. 
The random variable 9\ is time of the first velocity alternation, # 2 is time between 
the first and the second velocity alternations and 9 3 is time between the second and 
the third velocity alternations. 
We have that P {i — | < 6\ < t\ = e~ t+ i — e _ * and as it easy to calculate 

/ P \o 3 > t - u, 9 2 <^,e l e du) = (l - e~i) [ e~ t+u e~ u du 
Jo ^ 2 J J Q 

= (l-e-i)te- t . 
Whence, it is easily verified that lim e j_ {o<t-x(t)<e} _ 1 ^ g -t _|_ ie -ty 
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Similarly, P {t + x (t) < e} = P {9 2 > t - §, 0i < §} + o (e) . This implies that 

lim PQ + ,(t)< £ } = l e _ t=lim/c( ^ ) _ 

eio s 2 x\,-t 

Therefore, f c (t, x) is a solution of the Goursat problem for the linear second 
order hyperbolic equation that ensures the uniqueness of the solution of Eq.(12) 
with conditions (4). It means that / (t,x) is the pdf of the particle's position for 
to = 1. 
It is relevant to remark that function f (£, x) coincides with the result obtained in 

[5] \t 



Now, we turn to the case to — 2 and continue to calculate integrals of u[ 



It follows from u° A = 2-gf — 2uq that 



t 
u° 4 dy =2 

-t 



( 7T / u \dy — u\ (t, t) — u\ [t, —t) I — 2sinh t — 2sin t 
4(sinhi +sint -£)-2sinh£ -2sin£ =2sinht +2sin£ - U. 



p, 3 

Next, it follows from u\ = 2-gf — 2uq that 

/ u\dy =2 ( — / uldy - u\ (£, t) - u\ (t, -t) ) - 2sinh t + 2sin t 
=4sinht — 4sint — 2sinh£ + 2sin£ = 2sinh£ — 2sint . 
For t < \y\ we introduce the function g (t, y) — g c (t, y) + g s (£, y), where 

9c (t, y) = -ul (t, y) + - (u\ (£, y) + «f (t, y) + u\ (£, y) + u\ (i, y) + u° 4 (t, y)) , 

g s (t,y) = 6(t-y) + t6(t-y). (13) 

By construction, the function g c (£, y) is a solution of the following equation 

WW) 3(^y)-9(t,y) = 0. (14) 

Therefore, the function f c (t,x) — e' t g c (t,x) is a solution of Eq.(l) for to = 2 
(X = v = l). 

We put / (t, x) = f c (t, x) + e t g s (t, x). Taking into account the values of inte- 
grals of functions, which are involved in the expression for g c (t,y), we have that 
f_ t f (t, x)dx = 1, for all t > 0. 

Let us prove that lim xtt f c (t, x) = lim £ ^ P{ ° < '~ :c(t)<£} and \im xi - t fc {t, x) = 
lim 4 „ MM. 
It follows from Lemma 2 that for to = 2 we have 



and 



,. P{0< t-x(t) <£} 1 _ t 
hm — - — = -te . 

eio e 2 



hm nt + xitXe] = a 



eiO e 

It easily verified that lim^t u° (t, y) = 0, lim^ Uq (t, y) = and consequently 
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■iv ■"'■"--" nv2y/w=tf h y^ v > 2' 



lim g c (t, y) = lim —== = h ( \A 2 - y ] = -■ 



lim g c (*, j,) = lim Zl^ /i (V* 2 " 2/ 2 ) = 0. (15) 



vi-t vi-t 2^/t 2 - y 



Thus, 



, / x 1 t ,- P{t-x(t) <e\ 

lim f c (t, x) = -te~ l = lim — i ^ J - , 

xt* 2 eio e 

lim / c (i, x) = = lim P{t + ,T ^ <£} . (16) 

Let us show that conditions (15) with the condition J g (t, y) e~ t dx — 1 insure the 
uniqueness of the solution g c (t,y) for Eq.(14) and consequently, the uniqueness of 
the solution f c (£, x) of Eq.(12). 

It is easily seen that each solution of Eq.(ll) is a solution of Eq.(14). By changing 
the variables s = t + y, p = t — y, we reduce Eq.(14) to 

J^G{8,p)-G{8,p)=0. (17) 

Passing to the Fourier transform G (s, a) = J^° G (s,p) e lap dp in Eq.(17), we get 
the ordinary differential equation of order 4. Taking into account that lim^-j g c (t, y) 
0, we have 

G(0,a) = 0. (18) 

Hence, at most four independent solutions of the ordinary differential equation 
satisfy the initial condition (18) for each a. Passing to the inverse Fourier transform, 
we have four independent solutions of Eq.(14) with the condition lim^-t g c (t, x) — 
and just two of them satisfy Eq.(14) but not Eq.(ll). By construction, one of 
these solutions g c (t,y) is given by Eq.(13). As another solution we can take 

92(t,y) = ul(t,y) + u° 4 (t,y). 

It is easily verified that no linear combination c(t, y) of functions g c (t, y) and g2(t, x) 
satisfies conditions (16) and J (c(t,x) + g s (t,y)) e^ t dx = 1 for all t > 0, but 
solution g c (t, y). 

Therefore, the function / (t, x) is the pdf of the particle position at time t for m = 2, 
v = A = 1 and has the following form 

V ' 2 2Vt 2 - x 2 \ 



x 2 e * 



(ii [vt 2 - x 2 ) + Jx (yw 



2^(t 2 - x 2 f 



+5(t-x)e l +tS(t- x)e" 
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In much the same way as the pdf / (i, x) of the particle position for m = 2 was 
obtained we can also get solutions of Eq.(l) with conditions (2) and (4) for each 
m > 2. 
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